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Employing first principles total energy calculations we have 
studied the behavior of Ga and Al adatoms on the GaAs(OOl)- 
132 surface. The adsorption site and two relevant diffusion 
channels are identified. The channels are characterized by 
different adatom-surface dimer interaction. Both affect in a 
novel way the adatom migration: in one channel the diffusing 
adatom jumps across the surface dimers and leaves the dimer 
bonds intact, in the other one the surface dimer bonds are 
broken. The two channels are taken into account to derive ef- 
fective adatom diffusion barriers. From the diffusion barriers 
we conclude a strong diffusion anisotropy for both Al and Ga 
adatoms with the direction of fastest diffusion parallel to the 
surface dimers. In agreement with experimental observations 
we find higher diffusion barriers for Al than for Ga. 

PACS numbers: 68.35. Fx, 68.35 Ja 

Growth techniques, such as molecular beam epitaxy, 
operate under conditions far away from thermodynamic 
equilibrium. Particularly for growth at low tempera- 
tures or for structures with length scales smaller than the 
adatom diffusion length, features driven by the growth 
kinetics have been observed. In the GaAs/AlAs het- 
eroepitaxy the differences between the growth kinetics of 
AlAs and GaAs have been utilized to create low dimen- 
sional structures. Q Recently Kapon et al. || success- 
fully fabricated quantum wire heterostructure in which 
stimulated emission has been observed. 

Despite these successes the underlying microscopic 
processes such as adsorption, surface diffusion, desorp- 
tion, and nucleation are poorly understood. A key mech- 
anism in growth is cation surface diffusion which is con- 
sidered to be a rate limiting process. QJ^] Experimentally 
the surface diffusion is difficult to access. The deduced 
migration barriers ]o]-|To| for Ga adatoms span a range 
between 1.1 and 4.0 eV. Even the anisotropy of surface 
diffusion on the (001) surface is controversially debated: 
Shitara et al. [jllj speculate that fastest diffusion is along 
the [110] direction; in contrast, Kawabe et al. Q pro- 
pose the [110] direction. We have therefore performed 
first-principles total-energy calculations which constitute 
a powerful tool to study surface diffusion. 

Adatoms on solid surfaces occupy well defined bind- 
ing sites. The migration of the chcmisorbed adatoms can 
be described as a hopping between these sites. The ac- 
tivation energies for the individual hops are determined 



by the energy differences between the binding and the 
transition sites. These positions can be identified as the 
minima and saddle points of the potential energy surface 
over the configurational space spanned by the the coor- 
dinates of the adatom and the substrate atoms. In order 
to find all minima and transition sites a mapping of the 
entire configurational space is in principle required. How- 
ever, this is computationally neither possible nor useful. 
For the study of surface diffusion the mapping is com- 
monly restricted to a subspace given by the lateral coor- 
dinates of the adsorbate. Jl^l6[] This mapping gives the 
potential energy surface (PES) E(x, y) for a given lateral 
position (x, y) of the adatom where all substrate atoms 
and the z coordinate of the adatom are fully relaxed. 
This implies that all minima and saddle points relevant 
for the description of surface diffusion can be found on 
one single PES. However, in this Letter the results of our 
first principles study on the cation diffusivity show that 
this assumption fails for the migration on the GaAs(OOl) 
surface. 

In this paper we focus on the GaAs(001)-(2x4) sur- 
face in the (32 phase. This structure has been shown 
to be stable at the equilibrium |l^] as well as to be the 
dominating surface structure over a wide range of growth 
conditions. |]l8| , |l9f As shown in Fig. |l|, each (2x4) unit 
cell consists of two As dimers and two missing dimers in 
the topmost layer, a missing Ga pair in the second layer 
and a As dimer in the third layer. 

Our calculations of the PESs employ density func- 
tional theory in the generalized gradient approxima- 
tion (GGA). |2(J The electron- ion interaction is de- 
scribed by fully separable, norm-conserving pseudopoten- 
tials. ]^l],^2) The adatom-substrate system is modeled by 
a supercell geometry with a (4x4) periodicity parallel to 
the surface. This supercell has been tested to be suffi- 
ciently large to have a negligible adatom-adatom interac- 
tion. Perpendicular to the surface our supercell contains 
a vacuum of 6 layers GaAs and an additional layer of 
pseudo H-atoms (Z=0.75) to saturate the bonds of the 
lower surface. p3|] To prevent a dipole-dipole interaction 
between the inequivalent upper and lower surface of our 
slab we use a dipole correction. [Q The wave functions 
are expanded in a plane-wave basis with a cutoff energy 
of 10 Ry. The k-space integration was performed with 
a special k-point set, with a density equivalent to 64 k- 
points in the Brillouin zone of the (lxl) surface cell. 
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The PES has been mapped on an equidistant grid with 
a spacing of 1 A along the [110] and the [TlO] direction. 
At each position the upper four substrate layers and the 
adatom height were fully relaxed, starting from the posi- 
tions of the clean surface for the substrate and an adatom 
height of 3 A above the surface. In a second step we in- 
terpolated the PES from the energy values on the mesh. 
Subsequently we repeated our mapping around the in- 
terpolated saddle points using a finer mesh. The exact 
positions of the local minima were obtained by putting 
the adatom on the interpolated minima and relaxing all 
atomic positions without constraints. 

The resulting PES for a Ga adatom is shown in Fig. 
0(a). In each (2x4) unit cell we find the binding sites: 
Ai and the two equivalent positions A 2 and A 2 /. All 
three binding sites are on long bridge positions between 
surface As-dimer. This agrees with previous first princi- 
ples calculations Jl5|] for the metastable (2x4)-/3 struc- 
ture (assuming near-equilibrium growth conditions [p5| , 
at T = 900 K, a typical growth temperature, the surface 
area covered with the (i2 phase and the area covered by 
the P phase differ by more than two orders of magnitude) . 
At the three binding sites our calculations yield a strong 
inward relaxation of the adatom which resides almost 
in the center between the neighboring As atoms. The 
substrate atoms display only minor deviations (< 0.2 A) 
from their clean surface positions; in particular the sur- 
face dimer bonds remain intact. 

From PES shown in Fig. 0(a) migration on GaAs(OOl) 
can be described as a sequence of random hops on a pe- 
riodic (2x4) lattice with the three binding sites. The 
effective diffusion barrier is determined by applying the 
formalism of continuous time random walk (CTRW) ]26| ] 
that yields the diffusion tensor as a function of the indi- 
vidual hopping rates between neighboring binding sites. 
Within transition state theory |27) the individual hop- 
ping rate from a site Ai to a site Aj via a transition 
state Tk can be written as 

Tij = r§ exp[(£(T k ) - E(Ai))/k B T] (1) 

where £(Aj) and £/(Tk) are the energies given in Tab. @. 
r?j is a temperature independent prefactor and assumed 
to be roughly the same for all individual hops. This as- 
sumption leads to only minor errors as long as solely the 
effective diffusion barrier are discussed. @Jl|Jl| Apply- 
ing the CTRW formalism |28j and considering the size of 
the individual activation energies we find for the diffusion 
constants along the [110] and [TlO] directions 

n, , - r 12 r 21 (2r 22 +r 22 ,) „ 2 _ r i2 s 2 (9 s 
^[110] - (2r 12 +r 21 )(4r 22 +r 21 +2r 22 ,) Sa o~ — » a o \ z ) 

n - — 4ri 2 r 22 +rnr 2 i 9 2 ^ p 9 2 fo\ 
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ao is the lattice constant of GaAs. From the approx- 
imated solutions we get as effective diffusion barriers 



AE [1W] = E(T 1 )-E(A 1 ) and AE [JlQ] = S(T 4 )-£;(A 1 ), 
i.e. AE [1W] = 0.8 eV and A% 10] = 0.6 eV. 

In the above discussion we have implicitly assumed 
that all relevant minima in the configuration space could 
be found by relaxing the adatom from a position in the 
vacuum and the substrate from its ideal position in the 
clean case. This is true if all other minima of the total 
energy are either energetically much higher or separated 
from the first minima by a large barrier. The assump- 
tion should be correct when the surface exhibits mainly 
bulk-like bonds, as it is the case for most metal surfaces. 
For reconstructed semiconductor surfaces, however, sur- 
face bonds exist which differ significantly from the bulk 
bonds, as for example anion-anion dimers. On those sur- 
faces stable adsorption sites for a cation could be realized 
by breaking these weak surface bonds and forming bulk- 
like cation-anion bonds instead. 

Since in the PES shown in Fig. 0(a) the adatom does 
not break any surface As-dimer, we investigate in more 
details the adatom-dimer interaction. Fig. 0(a) shows 
the binding energy of the adatom as a function of its 
height z a -d on the dimer center with its lateral coor- 
dinates (x, y) fixed above the dimer center [site A3 in 
Fig. 0(b)] and the substrate fully relaxed. Approaching 
the surface from the vacuum we find a first local mini- 
mum at z a _d ~ 2 A which corresponds to a saddle point 
in the PES shown in Fig. 0(a). Coming closer to the 
surface there is a barrier at z a ^d ~ 1-2 A and another 
much deeper minimum at z a -d ~ 0.1 A. Such a double 
potential well was already predicted for Si on Si (001) 
by Roland and Gilmer. [|l6| The origin of the barrier 
and the presence of the second minimum can be under- 
stood by analyzing the bonding situation of the adatom 
to the dimer. At the first minimum the adatom interacts 
with the completely filled dangling bonds of the dimer 
[Fig. 0(c)], at the second one it breaks the dimer bond 
and forms directional bonds with the adjacent anions 
[Fig. 0(b)]. The barrier between the minima is due to 
the energy cost required to break the As-dimer. The 
binding energy at the second minimum is higher than 
that at the previously found binding sites Ai, A 2 and 
A 2 /. This shows that the adsorption site which is de- 
fined as the site with the highest binding energy cannot 
be on the PES displayed in Fig. |l](a). 

To find the real adsorption site and to include the in- 
fluence of the dimer breaking on the diffusion proper- 
ties we map the PES a second time. In contrast to the 
mapping of the first PES the adatom height z is now 
relaxed starting from a position 0.5 A above the surface 
and the surface As atoms are relaxed from a initial po- 
sition where the As- As dimer bond length is 3.5 A, i.e. 
the directional dimer bonds are broken. Fig. 0(b) dis- 
plays the resulting PES. Figs. |Tj(a) and (b) show that 
the two PESs are similar over a wide area but they dif- 
fer significantly around the dimers. In addition to the 
binding sites already present on the first PES we find 
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three new binding sites in each unit cell: A3 and the two 
equivalent positions A4 and A4/ . The comparison of the 
binding energies of the whole set of binding sites yields 
as adsorption site A3 , a short bridge site between the As 
atoms of one dimcr. 

At A3, A4 and A4/ the adatom breaks the dimcr, and 
charge accumulations between the adatom and the adja- 
cent anions confirm the formation of directional bonds. 
The angle between the two adatom-As bonds is about 
175°, i.e. the adatom is almost incorporated into the As 
top layer. This is accompanied by a strong horizontal 
relaxation of the As atoms bonded to the adatom. They 
relax along the [TlO] direction, and the As-As distance 
increases from w 2.5 A (the As-dimer bond length) to 
4.4 A which is close to the As-As distance in bulk GaAs 
(4.0 A). The relaxation keeps the nearest neighbor bond 
length of the substrate atoms nearly unchanged and re- 
sults in bulk like adatom-As bond lengths. Such a relax- 
ation has been observed experimentally for an Al termi- 
nated (2x1) reconstructed GaAs (001) surface. Q The 
driving mechanism for the flat incorporation of the cation 
is analogous to that which leads to the well known inward 
relaxation of the surface Ga atoms on GaAs(llO). The 
two-fold coordinated adatom with only partially occu- 
pied dangling bonds rehybridizes to a planar sp or sp 2 
bonding configuration. At the same time the three-fold 
coordinated As atoms bonded to the adatom gain energy 
from rehybridization to a p 3 bonding configuration. 

To find the exact binding energies at the transition 
sites T5, Tg and T7 between the short bridge and the 
long bridge binding sites the PES has been mapped 
around the interpolated transitions sites as a function 
not only of the lateral adatom coordinates but also of 
the dimer atoms. Note that we found an energy bar- 
rier to break the dimer for every pathway of the adatom 
toward A3 , A4 , and A4' . 

Including the new binding sites, surface diffusion is de- 
scribed by a random hopping on a periodic (2 x 4) lattice 
with six sites in each unit cell: Ai, A2, A2', A3, A4, 
and A4'. With the CTRW formalism we approximate 
the diffusion constants (the rather long exact expression 
will be published elsewhere (28|) as: 

T12T31 + Ti3r32 2 
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From theses equations we find as effective diffusion bar- 
riers AE [no] = £(T0 - E(A 3 ) and A£ [Tl0] = £(T 4 ) - 
E(Az). Inserting the energies given in Tab. [j| we find 
for the Ga adatom diffusion A£ , [ 110 ] = 1.5 eV and 
1.2 cV. From these values a significant dif- 
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fusion anisotropy with fastest diffusion parallel to the 
surface dimers can be concluded. 

The inclusion of the binding sites of the second PES 
changes the diffusion barriers dramatically. At first 



glance, taking into account only the second PES seems to 
be sufficient for the calculation of the diffusion barriers. 
However, from Eq. (||) it is clear that the description of 
the diffusion parallel to the surface dimers demands hop- 
ping rates from both PESs. If the contributions from one 
PES are neglected [as in the Eqs. (||) and (||)] , wrong bar- 
riers are obtained. The reason for involving both PESs 
is that the energy barrier to hop from a site Ai to a 
site A3 (which requires the breaking of a dimer) is only 
slightly higher than the one to migrate from the site Ai 
to a neighboring Ai site. Thus, once an adatom occupies 
a weakly bound Ai site it performs a number of hops to 
neighboring Ai sites (like surfing on the first PES) before 
being trapped again in a strongly bound A3 site. 

Our calculated diffusion barriers for a Ga adatom are 
at the lower limit of the experimentally deduced diffusion 
barriers (1.1 - 4.0 eV). This seems reasonable because the 
experimental results were derived rather indirectly from 
growth experiments and therefore they are affected by 
adatom-adatom and adatom-step interactions. Further, 
the adatom mobility depends on the As supply and de- 
creases with increasing As flux. ]30| ] Accounting for such 
effects should result in an effective diffusion barrier higher 
than that calculated here for a single adatom. 

We have performed similar calculations also for an Al 
adatom. The PESs are almost identical to those of a 
Ga adatom, i.e. the positions of the minima and of 
the saddle points for the two cations differ by less than 
0.1 A. The main difference is a slightly stronger corru- 
gation for Al adatoms leading to higher diffusion barri- 
ers: A_E[ 110 ] = 1.6 eV and Ai?p- 10 ] = 1.3. The enhanced 
corrugation can be understood in terms of the higher co- 
hesive energy of AlAs (£a°a s = 3.78 eV) compared to 
GaAs (^G a h As = 3.26 eV) implying that Al-As bonds are 
stronger than Ga-As bonds. From the higher barriers we 
conclude a lower mobility of Al adatoms which is in ac- 
cordance with several experiments. jlJU Our calculated 
difference in the barrier height of about 0.1 eV for Al com- 
pared to Ga adatoms agrees rather well with the value of 
0.16 eV obtained by Shitara et al. ^ from Monte-Carlo 
simulations of RHEED measurements of the growth of 
AlAs and GaAs. 

In conclusion, we have shown that the adatom-surface 
dimer interaction is crucial for determining the adsorp- 
tion site as well as the effective diffusion barrier. The ori- 
gin of the very stable new adsorption site is that a weak 
surface bond (As-As) is replaced by almost bulk-like 
cation-As bonds. This mechanism should also work on 
other semiconductor surfaces exhibiting surface dimers 
or trimers. 

Stimulating discussions with E. Pehlke and K. Shi- 
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nancial support. 
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TABLE I. Binding energies (eV) of a Ga adatom at various 
surface sites on the GaAs(001)-(2x4) /32 structure 
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FIG. 1. Total energy surfaces for a Ga adatom on the 
GaAs(001)-(2x4) (32 surface, (a) PES obtained when the 
adatom is relaxed from 3 A above the surface (b) PES ob- 
tained when the adatom is relaxed from 0.5 A above the sur- 
face with the surfac e dimers initially broken. The dashed box 
shows the unit cell. The contour-line spacing is 0.2 eV. The 
atomic positions of the clean surface are indicated for atoms 
of the upper two layers and for the As-dimers in the third 
layer (As: empty circles, Ga: filled circles). 
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FIG. 2. (a) Binding energy of an Ga adatom as function 
of its height z a -d above the center of an As-dimer. (b) and 
(c) charge density cross section through the adato m and the 
dimer atoms at the two minima of the binding energy. 
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